function M = getMomentMatrix(lattice)

switch lattice
    
    case 'D2Q9'
   
        M = [1 1 1 1 1 1 1 1 1;
            -4 -1 -1 -1 -1 2 2 2 2;
            4 -2 -2 -2 -2 1 1 1 1;
            0 1 0 -1 0 1 -1 -1 1;
            0 -2 0 2 0 1 -1 -1 1;
            0 0 1 0 -1 1 1 -1 -1;
            0 0 -2 0 2 1 1 -1 -1;
            0 1 -1 1 -1 0  0 0 0;
            0 0 0 0 0 1 -1 1 -1];
        
    case 'D3Q15'
        M=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1;
            -2 -1 -1 -1 -1 -1 -1 1 1 1 1 1 1 1 1;
            16 -4 -4 -4 -4 -4 -4 1 1 1 1 1 1 1 1;
            0 1 -1 0 0 0 0 1 -1 1 -1 1 -1 1 -1;
            0 -4 4 0 0 0 0 1 -1 1 -1 1 -1 1 -1;
            0 0 0 1 -1 0 0 1 1 -1 -1 1 1 -1 -1;
            0 0 0 -4 4 0 0 1 1 -1 -1 1 1 -1 -1;
            0 0 0 0 0 1 -1 1 1 1 1 -1 -1 -1 -1;
            0 0 0 0 0 -4 4 1 1 1 1 -1 -1 -1 -1;
            0 2 2 -1 -1 -1 -1 0 0 0 0 0 0 0 0;
            0 0 0 1 1 -1 -1 0 0 0 0 0 0 0 0;
            0 0 0 0 0 0 0 1 -1 -1 1 1 -1 -1 1;
            0 0 0 0 0 0 0 1 1 -1 -1 -1 -1 1 1;
            0 0 0 0 0 0 0 1 -1 1 -1 -1 1 -1 1;
            0 0 0 0 0 0 0 1 -1 -1 1 -1 1 1 -1;];
        
    case 'D3Q19'
        M=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1;
            -30 -11 -11 -11 -11 -11 -11 8 8 8 8 8 8 8 8 8 8 8 8;
            12 -4 -4 -4 -4 -4 -4 1 1 1 1 1 1 1 1 1 1 1 1;
            0 1 -1 0 0 0 0 1 -1 1 -1 1 -1 1 -1 0 0 0 0;
            0 -4 4 0 0 0 0 1 -1 1 -1 1 -1 1 -1 0 0 0 0;
            0 0 0 1 -1 0 0 1 1 -1 -1 0 0 0 0 1 -1 1 -1;
            0 0 0 -4 4 0 0 1 1 -1 -1 0 0 0 0 1 -1 1 -1;
            0 0 0 0 0 1 -1 0 0 0 0 1 1 -1 -1 1 1 -1 -1;
            0 0 0 0 0 -4 4 0 0 0 0 1 1 -1 -1 1 1 -1 -1;
            0 2 2 -1 -1 -1 -1 1 1 1 1 1 1 1 1 -2 -2 -2 -2;
            0 -4 -4 2 2 2 2 1 1 1 1 1 1 1 1 -2 -2 -2 -2;
            0 0 0 1 1 -1 -1 1 1 1 1 -1 -1 -1 -1 0 0 0 0;
            0 0 0 -2 -2 2 2 1 1 1 1 -1 -1 -1 -1 0 0 0 0;
            0 0 0 0 0 0 0 1 -1 -1 1 0 0 0 0 0 0 0 0;
            0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 -1 1;
            0 0 0 0 0 0 0 0 0 0 0 1 -1 -1 1 0 0 0 0;
            0 0 0 0 0 0 0 1 -1 1 -1 -1 1 -1 1 0 0 0 0;
            0 0 0 0 0 0 0 -1 -1 1 1 0 0 0 0 1 -1 1 -1;
            0 0 0 0 0 0 0 0 0 0 0 1 1 -1 -1 -1 -1 1 1;];
        
        
        
    otherwise
        disp('Unknown lattice type');
       
        
        
end